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RANDOM ORIENTED TREES: A MODEL OF DRAINAGE 

NETWORKS 

By Sreela Gangopadhyay, Rahul Roy^ and Anish Sarkar^ 
Indian Statistical Institute 

Consider the d- dimensional lattice Z'^ where each vertex is "open" 
or "closed" with probability p or I — p, respectively. An open vertex 
V is connected by an edge to the closest open vertex w such that 
the dth co-ordinates of v and w satisfy w{d) = v{d) — 1. In case of 
nonuniqueness of such a vertex w, we choose any one of the closest 
vertices with equal probability and independently of the other ran- 
dom mechanisms. It is shown that this random graph is a tree almost 
surely for d—2 and 3 and it is an infinite collection of distinct trees 
for d > 4. In addition, for any dimension, we show that there is no 
bi-infinite path in the tree and we also obtain central limit theorems 
of (a) the number of vertices of a fixed degree v and (b) the number 
of edges of a fixed length 

1. Introduction. Leopold and Langbein (1962) introduced a geometric 
model of natural drainage network which they described as 

using a sheet of rectangular cross-section graph paper, each square is presumed 
to represent a unit area. Each square is to be drained, but the drainage channel 
from each square has equal chance of leading off in any of the four cardinal 
directions, subject only to the condition that, having made a choice, flow in the 
reverse direction is not possible. Under these conditions it is possible for one 
or more streams to flow into a unit area, but only one can flow out. 

Subsequently Scheidegger (1967) introduced a direction of flow. In his study 
of Alpine valleys, he imposed conditions on the Leopold and Langbein model 
by requiring that the drainage paths be in the "direction of high gradients 
between watershed and main valleys." Thus the drainage forms an oriented 
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network, with a square emptying to one of its two neighbors in a preferred di- 
rection. Howard (1971) removed the restriction of drainage to a neighboring 
square and modehed a network to include "headward growth and branching 
in a random fashion." Rodriguez-Iturbe and Rinaldo (1997) present a survey 
of the development of this field. 

The random graph we study here follows the one described by Howard 
(1971) with the caveat that a stream is not permitted to terminate or become 
inactive. Thus we consider the d-dimensional lattice %^ where each vertex is 
"open" or "closed" with probability p or 1— p, respectively. The open vertices 
represent the water sources. An open vertex v is connected by an edge to 
the closest open vertex w such that the dth co-ordinates of v and w satisfy 
w[d) = v{(£) — \. In case of nonuniqueness of such a vertex we choose any 
one of the closest vertices with equal probability and independently of the 
other random mechanisms. These edges represent the channels of flow in the 
drainage network. 

Our main result (Theorem 2.1) is that, for c? = 2 and 3, all the tributaries 
connect to form one single delta, while for d > 4, there are infinitely many 
deltas, each with its own distinct set of tributaries. In this connection it 
is worth noting that (Theorem 2.2) there is no main river, in the sense 
that there is no bi-infinite river; instead, each tributary has its own distinct 
source. In addition, for any dimension, we obtain central limit theorems of 
(a) the number of sites where a fixed number v of tributaries drain, as well 
as of (b) the number of channels of a fixed length I. 

Similar tree-forest dichotomies have been studied for the uniform span- 
ning tree model by Pemantle (1991) and for the minimal spanning tree model 
by Newman and Stein (1996). Ferrari, Landim and Thorisson (2002) have 
obtained similar results for a continuous version of this model. 

In two dimensions we obtain the main result by showing that the distance 
between two streams starting at two different sites forms a martingale and 
thereby invoking the martingale convergence theorem. For three dimensions 
we employ a technique based on Lyapunov functions, while in four or higher 
dimensions we couple the streams starting at two different sites with two 
independent and identically distributed random walks starting at these two 
sites. To show that there are no bi-infinite paths in the graph we utilize 
the stationarity of the model and use a Burton-Keane type argument. The 
limit theorems are obtained by checking that the random processes satisfy 
the conditions needed to apply Lyapunov's central limit theorem. 

The formal details of the model and the statements of results are in the 
next section. 

2. The model and statement of results. Let = {0, 1}^^ and let T be 
the a algebra generated by finite-dimensional cylinder sets. On (0,jr) we 
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assign a product probability measure Pp which is defined by its marginals 
as 

Pp{uJ : uj{u) = 1} = 1 — Pp{uJ : uj{u) = 0} =p 

for n G Z"' and < p < 1 . 

Let {Uu,v - UjV & Z'^,v{d) = u{d) — 1} be i.i.d. uniform (0, 1] random vari- 
ables on some probability space Here and subsequently we express 
the co-ordinates of a vector u as u = (ti(l), • ■ • ,u{d)). 

Consider the product space x E, x S ,¥ := Pp x fi). For {io,C) & ^xE, 
let V(= V(u;,^)) be the random vertex set defined by 

V{uj,0 = {ueZ'^:uj{u) = 1}. 

Note that if u G V(ti;,^) for some G H, then u G V(ci;,^') for all ^' G H and 
thus we say that a vertex u is open in a configuration u; if n G V(a;,^) for 
some ^ G H. 
For u G Z"', let 

Nu = Nu{io,C) 




V G V{uj, : v{d) = u{d) - 1 and 



d ( d 

\ v{i) — u{i)\ = min< ^ \ w{i) — u{i) \ : w G V(w, 
i=l U=i 

w{d) = u{d) - ijj 

Note that for p> 0, is nonempty almost surely and that is defined 
for all u, irrespective of it being open or closed. For u G Z'^, let 

h{u) G Nu{lo,^) be such that 

Again note that for p> and for each u G Z'^, h{u) is open, almost surely 
unique and h{u){d) = u{d) — 1. On V(u;,^) we assign the edge set 8 = 

£{u^,0 :={(n,/i(u)):nGV(u;,0}- 

Consider that graph Q = (V,<f) consisting of the vertex set V and edge 
set £. For p = 0, V = almost surely, and, for p = 1, (u, u ) G <5 if and only if 
u{i) = v{i) for alH / d and \u{d) — v{d)\ = 1. Also, for a vertex u G V(ti;,^), 
there is exactly one edge "going down" from u; that is, there is a unique edge 
{u,v) with v{d) < u(d). Thus the graph Q contains no loops almost surely. 
Hence, for <p < 1, the graph Q consists of only trees. Our first result is 
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Theorem 2.1. Let <p < 1. For d = 2 and d = 3, G consists of one sin- 
gle tree W-almost surely, while for d> 4, Q is a forest consisting of infinitely 
many disjoint trees ¥-almost surely. 

Regarding the geometric structure of the graph Q, we have 

Theorem 2.2. Let < p < 1. For any d>2, the graph Q contains no 
bi-infinite path ¥-almost surely. 

Now for v >0, let Sn be the number of vertices in V H ([l,?!]"^) of the 
graph Q with degree v +1. Also, for I > 1, let L„ be the number of edges of 
Li-length I in the graph Q with one end vertex in V H ([l,n]'^). 

Theorem 2.3. Asn^oo: 
i^) '^""^fi^ converges weakly to a normal random variable; 
(b) '^""^2^""* converges weakly to a normal random variable. 

Finally, for d = 2, given that a vertex v is open, the following proposition 
gives the exact distribution of the degree of v. 

Proposition 2.1. Given that a vertex v is open, the degree of the vertex 
in the graph Q has the same distribution as that of 1 + Y -\- Xi + X2, where 
Y , Xi and X2 are independent nonnegative random variables such that 



Thus the expected degree of a vertex, given that it is open, is 2. 

Remark 2.1. As in Lemma 7 of Aldous and Steele (1992), using the er- 
godicity of the process, it may be shown that in any dimension, the expected 
degree of a vertex, given that it is open, is 2. 

3. Proof of Theorem 2.1. We fix < p < 1 and for u,v E Z'^~^ consider 
the d-dimensional vectors u := (u, 0) and v := [v, 0) and let {X^, —n) := /i"(u), 
where h"' denotes the n-fold composition of h defined in (1). For Zn{= 
Zn{u, v)) := X^ — X", we first observe that it is a time-homogeneous Markov 
chain with state space Z'^"^; indeed, this follows on writing {Zn+i = Zn+i, 

Zn = Zn, . ■ . , Zq = Zq} = [jx„+i,...,xoe'E'i{-^x^^ = ^n+l,X'^^^zo = + Zn+1, 
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with probability 1 — p 
with probability p, 




for r = 0, 



for r > 1. 
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= a;„, =Xn + Zn,...,X^^= xo, =Xo + zq} and using the 
Markovian property of the process -.n > 0}. 

The connectedness or otherwise of the graph Q is equivalent to whether 
or not Zn is absorbed at the origin. For d = 2 and 3, we show that Zn 
gets absorbed at the origin, G Z"^^^ with probabihty 1; while for d> 4, 
Zn is a transient Markov chain and hence has a positive probability of not 
being absorbed. In this connection observe that instead of the above Z„, 
if we had considered a modified Markov chain Z„, where is no longer an 
absorbing state, but from we move in one step to some fixed vertex u 7^ 
with probability 1 and the other transition probabilities are kept unchanged, 
then to show that the original process Zn is absorbed at almost surely, it 
suffices to show that the modified Markov process Zn is recurrent. A more 
formal argument for this would require Zn and Zn to be coupled together 
until they hit the origin, which occurs almost surely if the modified process 
is recurrent. For the case d = 3, we will show that Zn is recurrent. The proof 
is divided into three sections according as d = 2, d = 3 and d > 4. 

3.1. d = 2. Fix i < j and observe that X" < for every n > 1, where 
X" and X" are as defined earlier. Thus the Markov chain Z„ := X" — X" 
with Zq = j — i has as its state space the set of all nonnegative integers. Since 
the marginal distributions of the increments of X" and X" are identical 
with finite means, {Zn : n > 0} is a nonnegative martingale. Hence, by the 
martingale convergence theorem [see Billingsley (1979), Theorem 35.4, page 
416], Zn converges almost surely as n — > 00. Since {Z„, : n > 0} is also a time- 
homogeneous Markov chain with as the only absorbing state, we must have 
Zn ^ as n ^ 00 with probability 1. Since this is true for all i < j, we have 
the result for d = 2. 

3.2. d = 3. Throughout this section the letters u, v in bold font denote 
vectors in Z'^, u, v in roman font denote vectors in and u, v in italic font 
denote integers. Fix two vectors u := (u, 0) and v := (v, 0) in 7? x {0} and let 
Zn{= Z„(u, v)) be the time-homogeneous Markov chain with state space Z^ 
as defined at the beginning of this section. We shall exhibit, by a Lyapunov 
function technique, that this Markov chain Zn is recurrent, thereby showing 
that Zn is absorbed at the origin with probability 1. 

Consider the function /:M^ — > [0,oo) defined by /(x) := y^log(l + HxHl) 
where || • II2 is the standard L2 norm (Euclidean distance). Since /(x) — > 00 
as ||x||2 00, by Foster's criterion [see Asmussen (1987), Proposition 5.3 of 
Chapter I, page 18] the following lemma implies that Z.„ is recurrent. 

Lemma 3.1. For all n>0, there exists T > such that, for all ||x||2 > T, 
we have 

E(/(Z„+i)-/(Z„)|Z„ = x)<0. 
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Proof. Let g : [0, oo) [0, cxd) be defined as g{x) := ^log(l + x). Clearly 
g{x) > for all X > and g{x) — > oo as x ^ cxo. Also, for x, y > 0, the Taylor 
series expansion yields 

(2) g{x) -g{y)<{x- y)g('^ (y) + ^-^^^g^'^ (y) + ^^^^^'^ iv) , 

which holds because the fourth derivative 
3 11 



^^'^^'^ {l + srg{s) 4(l + s)4(5(s))3 
18 15 

< for s > 0. 



8{i +sng{s)r ni+sng{s)y 

The first three derivatives of g, which we will be needing shortly, are 

1 



2(l + s)<7(s)' 
1 
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2{l + s)^gis) 4(l + s)2(5(s))3^ 



il + s^gis) 4(l + s)3(5(s))3 8(l + s)3(5(,))5- 
Note that, for all s large, 

Assuming for the moment that (we will prove this shortly), for some a > 0, 

(3) m\Zn+l\\l - WZnWl I = X) = a + o(||x||2 2), 

(4) - WZnWlf I = x) > 2a||x||^, 

(5) - WZnWlf I = X) = 0(||x||2), 

as ||x||2 oo, and using the above estimates and expression for derivatives, 
we have, for all [5 := ||x||2 large and for some nonnegative constants Ci and 
C2, 



E(/(Z„+i)-/(Z„)|Z„ = x) 

a + Ci//3 2a/3 



< 



2(l + /3)Vlog(l + /3) 4(l+/?)Vlog(l + /3) 
2a/3 ^ 3C2/3 



8(l + /3)V(log(l + /3))^ (l + /3)Vlog(l+/3) 
1 



4Ci 24C2/3 2a/3 
4a + 4Ci + — + 



8(l + /?)Vlog(l+/3)L /3 1 + /? log(l + /3) 
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The term inside the square brackets tends to — oo as /? — > oo; therefore, for 
all sufficiently large (3, the term is negative. Thus to complete the proof of 
the lemma we need to show (3)~(5). 

Let Dk '■= {v S : ||v||i < k} denote "Li-diamond" of radius k and 6Dk '■= 
{v G : ||v||i = A;} its boundary, where || • ||i denotes the Li norm. Consider 
the probability distribution of the step size of the random walk, associated 
with the tree generated by one particle, that is, the distribution of : 

p,:=P(Xi=u) 

(6) f p, ifu = o, 

(1 _ _ (1 _ p)#5£'fc) 

^ ^—^ iornG6Dk,k>l, 

where o := (0, 0) is the origin and i^A denotes the cardinality of the set A. 
For any k>l and i,j > 0, define 

mi{k):= u\pu 

u: = {ui,U2)eDk 

and 

mij{k):= u\ul,p^. 

u: = (iii,U2)eDfe 

Since {—ui,—U2) G whenever (ni,M2) G D^, it is clear that, for every 
A; > 1 we have 

, , rriiik) = for all odd i and 

(7) 

mij[k) = whenever either i or j is odd. 

Further, since #Dk = 1 + 2k{k + 1) and #5Dk = ik, we have that, for all 
even i, 

< mj :=limmj(A;) = ^ u\pu 

u:(ui,U2)eZ2 

< ^(#5Dfc)(max{^z, : (^,^2) G D^})'^- ^ j—^ ^ 

k=i *^^^ 

00 

= Y '^'(l - (1 -p)"') 

fc=i 

< 00. 

Similarly, when both i and j are even, mij{k) — > k 00, where < 

•= Z]ugz2 '^\'U2P\x < 00. Further, p^ being the same for every u on 6Dk, 
the various quantities nii and mj j remain unchanged if in their definitions 
we had considered U2 instead of ui. 
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Moreover, 

k'^ {m2 - in2{k)) < ^ ufpu 

oo 
j=k+l 

as /c ^ OO since the sum 

oo 

^/(1_p)1+2j0-i)<oo. 
A similar result holds for m,o(A;) and so we have 

(8) m2{k) = 1712 + o{k^'^) and mo(A;) = mp + o(A;^^) as/c— >oo. 
Now we proceed to compute the expectations: 

E(||Z„+i||2-||Z„||2|Z„=x) 

= E (||x + a-bi-||x||i) 

a,beZ2 

X E [IP{^u+' = K + a, 

WGZ2 

(9) X:+' = x: + h\X: = w,X: = w + x}] 

X [P{X; = w,X;^ = w + x|Z„ = x}] 

= E (||x + a-b||i-||x||2) 

X F{X^ = x + a,Xl=h\X^ = o,X^ = x}, 

where we have used the translation invariance of the model. 

To calculate the above sum we let k := ||x||2/4. Note, for a,b € 1?^, we 
have P{X^ = X + a, = b I X° = o, ^0 = x}0 = PaPb; thus, using (7) and 
(8), 

ri(l):= (IN + a-bi-||x||2) 



X P{X^ = x + a,Xo^ = b I X° = o,X° = x} 



(10) 
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= [(ai-^i)' + 2xi(ai-6i) 

+ (02 - ^2)^ + 2X2(02 - b2)]Pa.Ph 

= 4'm2{k)mo{k) 

= 4m2 + o(fc~^) as /c — > 00. 

Also, if b ^ Dk, then, taking ||b||i = k + 1 for some / > 1, the occurrence of 
the event {X^ = b} requires that ah the vertices in the diamond Dk+i-i be 
closed and that at least one vertex of 5D]^j^i be open — an event which occurs, 
with probability (1 - p)i+2{fc+«-i){fc+0 _ (1 _ p)i+2{fc+0{fc+«+i). Moreover, if 
{Xq = b} occurs, then X^^ must lie in the smallest diamond centered at x 
which contains the vertex b; thus \\X^ — X0II2 < H-'^xlli + ll-^olli ^ (ll^lli + 
l|b||i) + ||t>||i = 6A; + 21. Now noting that there are 4(A; + /) vertices on SD^j^i 
and that an argument similar to the above may be given when a ^ -D^, we 
have 

T2(l):= (I|x + a-b||i-||x||i) 

a^Dfe or b^Dfc 

X ¥{Xl = X + a, = b I X° = o,X° = x} 

(11) 

< 2 ^ 4(A: + /) ((6A: + 2lf + (4A;)2) (1 - p)i+2(fc+;-i)(fc+0 
X [1- (l-p)^('=+')] =o(A:-2) ask ^00. 



This establishes (3) with a = 4m2. 

For (4), calculations as in (9) show that E((||Z„+i||^ - \\Zn\\lf | Z„ = x) > 
Ti(2) where, performing calculations as in (10), 

Ti(2):= (l|x + a-b||2-||x||2)2p{xi=x + a,Xi=b|XO = o,xO = x} 

a,bGDfe 

= H [(ai -6i)2 + 2xi(ai -61) + (02 -62)^ + 2x2(02- 62)] Vpb 

= 8m2(A;)mo(A;)||x||2 + Ami{k)mQ{k) + 16(?Ti2(A;)) + 4m2,2(^)"io(^)- 

Now, as k ^ 00, from (8) we have m2(fc)mo(fc)||x||2 — m2||x||| = 16(?7i2 + 
o(fc-2))(l + o(A;-2))A;2-16m2A:2 = IQk"^ o{k-'^) + o{k-'^) = o(l); also 4m4(A;)mo(it) + 
16(m2(A;))^+4m2,2(fc)mo(A;) ^ 4m4 + 16m| + 4m2,2 > 0. This establishes (4). 
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Finally, for (5), we write E((||Z„+i||2 - \\Zn\\lY \ Z„ = x) = ri(3) + r2(3), 
where 

ri(3):= (l|x + a-b||l-||x||imxi=x + a,Xi=b|XO = 0,XO=x} 

= [(ai - fci)^ + 2xi(ai- 61) + (02-62)2 + 2x2(02- 62)] VPb 

= P^Ph[{ai-hif{2xi+ai-hif 

+ 3(01 - &i)2(2j;i + ai - bif{2x2 + 02 - 62)(a2 - ^2) 

+ 3(ai - 6i)(2xi + ai - 6i)(2x2 + 02 - 62)^(02 - 62)^ 

+ (a2-?'2)^(2x2 + a2-62)^] 
= ri,i(3) +ri,2(3) +ri,3(3) + ri,4(3) (say). 

Now, 

^1,1(3):= E PaPb(«i-6i)^[(ai-^>i)^ + 6x1(01-61)2 

a.beDfc 

+ 12x?(ai-6i) + 82;?] 

= X] PaPb(ai - 61)*^ + 12Xi E PaPb(oi-6l)^ 
a,b£Z)j; a,bsDfc 

= C3(A:) + C4(A:)x?, 

where C2,{k) and C4(/c) are both polynomials in mi(k) and mij{k), each of 
which converges to the corresponding polynomial in rrii and A; — > 00. 

Similar calculations show that Ti^4(3) := X)a,be-Dfc PaPb(«2 — 62)'^ [(^2 — 62)'^ + 
6x2(02 - 62)2 + 12xi(a2 - 62) + 8xi] = C-sik) + C4ik)xl 
Also, 

ri,2(3)=3 Yl P^Ph{ai-bif[{ai-hf + Axi{ai-bi) + 4xl] 

a,bGDj. 

X [2x2(02 - 62) + (a2 - 62)^] 

= 3 E PaPb(oi - 61)^(02 - 62)^ 

a,beDfe 

+ 12X? Y PaPb(«l- 61)^(02-62)^ 
a,b6Dfc 
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= C5{k)+Ceik)xl 

where, as above, C^{k) and CQ{k) are both polynomials in mi(k) and niij^k), 
each of which converges to the corresponding polynomial in rrii and rrii j 
as k ^ oo. Similar calculations show that Ti^3(3) := 3X]a bs-Di, ^'aPb(a2 — 
b2f[ia2 - h2? + 4x2(a2 - 62) + 4xi][2xi(ai -hi) + (ai - bif] = C^{k) + 
C^{k)xl 

Finally, calculations similar to (11) yield T2(3) := J^nifDk or b^Dfe(ll^ + ^~ 
b||i - ||x||i)3p{Xi =x + a,Xl = h\ = 0,X^=x} = o{k-^) as A; ^ 00. 
This establishes (5) and completes the proof of Lemma 3.1. □ 



3.3. d>4. For notational simplicity we present the proof only for d = 4. 
Throughout this section the letters u, v in bold font denote vectors in Z^, u, 
V in roman font denote vectors in and u, v in italic font denote integers. 
We first show that on Z^, the graph Q admits two distinct trees with positive 
probability, that is, 

(12) ¥{g is disconnected} > 0. 

Consider a random vector X G Z^ defined as follows: for /c > 0, let := 
{v G Z^ : ||v||i < /c} denote the three-dimensional diamond of radius k and let 
(5Afc := {v G Z'^ : ||v||i = A:} denote its boundary. As in (6), the distribution 
of the random vector X is given by 



p, if V 

(13) ¥{X = v) = <( (l-p)#^fc-i(l-(l-p)#^^fe) 



for V G 5Afc, > 1, 



where o := (0,0,0) and ij^A denotes the cardinality of the set A. It may 
easily be checked that X^veza P{X = v) = 1. 

Next, for a fixed vector u := (n(l), . . . , m(4)) G Z^, consider the graph Ti := 
(VU{u},£:u{(u,/i(u))}). For n > 0, let /i'^(u) := (ff"(u),t) for ff"(u) G Z^ 
and t = n(4) — n G Z. Here we take h^{\i) = u. Observe that for fixed u, g^{u) 
has the same distribution as (u(l), m(2), ii(3)) + J22=i-^i^ where Xi,X2,... 
are i.i.d. copies of X. Hence {g'"(u):n > 0} is a symmetric random walk 
starting at ^'^(u) = (u(l), m(2), u(3)), with i.i.d. steps, each step size having 
distribution X. However, for v G Z^ with v{4:) = m(4), in the graph (V U 
{u, v},£: U {(u,/i(u)), (v,/i(v))}) the processes {5"'(u)}„>o and {5"(v)}„>o 
are not independent and so, to obtain our theorem, we cannot use the fact 
that, with positive probability, two independent random walks on Z'^ do not 
intersect. Nonetheless, if u and v are sufficiently far apart, their dependence 
on each other is weak. In the remainder of this section we formalize this 
notion of weak dependence by coupling two independent random walks and 
the processes {<^"(u), <^"(v) :77, > 0} and obtain the desired result. 
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For V = (v, 0), given e > define the event 

(14) 

g\v)^g\0) foralli=l,...,n^}, 

where := (0, 0, 0, 0). 

Lemma 3.2. For < e < 1/3, there exist constants C,f5 >0 and no > 1 
such that, for all n>nQ, 

inf P(A„,,(v))>l-Cn-^. 
sO(v)eA„i+AA„i-e 

Assuming the above lemma, we proceed to complete the proof of (12). 
We shall return to the proof of the lemma later. 

For i>l and n > no, let Ti(= rj(n)) := 1 + n^ + (n^)^ H h (n"')^' and 

take To = 1. For fixed v, we define 

Bo = i?o(v) := {g(v) G g{0) + (A„i+. \ A„i-.)}, 

and having defined -Bo, ... , Bi^i, we define 

Bi = B,{^r) := {^-'(v) G 9^'{0) + (A^,.(,+,) \ A^,.,,_^,) 

and ^■'■(v) / ^■'■(O) for ah n^i + 1 < j < n}. 



Clearly, 



P{£,J(v) / g^{0) for all j > 1} 



> 



(15) 



.1=0 



limp n^?. 



i— >oo ■ , 

1=1 \ j=0 



Since ¥{Bo) > 0, from (15) we have that F{g^^r) / ^■'■(O) for all j > 1) > 
if E^,l-F{B,\n'-JoBj)<oo. 

For fixed I > 1, let Ui := /i^'(O) and vi := /i^'(v). Now {(/i"(0), /i"(v)) : n > 
0} being a Markov process and, since g^{-vi){uj,S^) G (7'^(ui)(u;, ^) + (A^2*(i+e) \ 
A„2i{i-e)) for (w,^) G -Bz(v), we have 
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> inf P{^^("')'' (vi) G (ui) + (A^,,+i(,+,) \ A^,,+i(,_,,), 

(16) /(ui) + g'^iy^) for all = 1, 2, . . . , (n^f } 

= infP(^„,,^(u)) >l-C(n2V^ 

where infi is the infimum over all Ui, vi G with ^^(vi) G g(''(ui) + (A^2'(i+e) \ 
^„2i{i-e)) and inf2 is the infimum over all u with 5° (u) G (A^jifi^,) \ A^2;(i_,)) 
and the last inequality follows from Lemma 3.2. Thus EiSi(l -I^C^/I n5=o ^i)) 
CZ]/^i(?^^ < 00, thereby completing the proof of (12). 

To prove Lemma 3.2, we have to compare the trees {/i"(0)} and {/i"'(v)} 
and independent "random walks" {0 + (X]"=i —n)} and {v + (X]"=i Yi, —n)}, 
where {Xi, X2, . ■ .} and {Yi,Y2, . . .} are independent collections of i.i.d. 
copies of the random variable X given in (13). 

We now describe a method to couple the trees and the independent ran- 
dom walks. Before embarking on the formal details of the coupling proce- 
dure, we present the main idea. 

From a vertex we construct the "path" {0 + {J^i^iXi,—n)}. Now 
consider the vertex v with v = {vi,V2,V3,0). In case the diamond D := 
{u G : ||u||i < ||Xi||i} is disjoint from the diamond D' := {u G : ||u — 
(ui,?;2,?;3)||i < ll^^lli}, then we take h^{v) = {v + (Yi, -1)}. If the two di- 
amonds are not disjoint, then we have to define /i^(v) taking into account 
the configuration inside the diamond D. Similarly, we may obtain /i^(v) by 
considering the diamonds {u G : ||u — Xi\\i < ||X2||i} and {u G Z'^ : ||u — 
5^(^)111 < ||i^2||i}- Note that if, for each i = 1, . . . ,n, the two diamonds in- 
volved at the ith stage are disjoint, then the growth of the tree {(/i*(0), /i*(v)) : 
i < n} is stochastically equivalent to that of the pair of independent "random 
walks" (0 + (ELi Xi, -n), V + (ELi Yi, -n)). 

We start with two vertices u := (u, 0) and v := (v, 0) in Z^ with u, v G Z'^. 
Let {C/f (z) : z G Z^},{U^{z) : z G Z^} and {U{{z) : z G Z^}, {U^{z) : z G Z^} be 
four independent collections of i.i.d. random variables, each of these random 
variables being uniformly distributed on [0, 1] . 

Let /cu and 1^ be defined as 

fcu := min{A; : C/{'(z) < p for some z G (u + A^)}, 
Iv ■= min{/ : C/^(z) < p for some z G (v + A/)}. 

Now define as 

mv := min{m : either Ui{z) < p for some z G (v + A^) \ (u + A^^^) 
or Ui{z) < p for some z G (v + A^) Pi (u + A^^)}. 
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Also, define the sets 



u 



{zG(u + AfeJ:[/i-(z)<p} 




{zG(v + A;J:C/i^(z)<p}, 

{z G (v + A^J \ (u + AfcJ : U^^) < p} 

U{zG (v + A^Jn(u + AfcJ:C/{'(z) <p}. 



We pick: 



(a) (j){u) G iVu sucli tliat U^{(p{u)) = mm{U^{7.) : z G iVu}; 

(b) C(v) G A^v^ sucli tliat C/2''(C(v)) = min{J72''(z) : z G N^}; 

(c) ^/;(v) G sucli that C/2''(V'(v)) = min{J7J(z) : z G A^}. 

Taking (/'^(u) = u, (j)"'{u) = (n)) , and similarly for C"(v) and ip^iy), 

we note that the distribution of {(((^"(u), — n), (<^"(v), — n)) : n > 0} is the 
same as that of {((u + X^iLi —n), (v + J2^=i ^i, —n)) : n > 0}, that is, two 
independent "random walks," one starting from (u,0) and the other start- 
ing from (v,0). Also the distribution of {(/i„(u, 0), /i„(v, 0)) :n > 0} and that 
of {(((/>" (u), — n), {ip"'{v), —n)) :n > 0} are identical. Thus, the procedure de- 
scribed above may be used to construct the trees from (u, 0) and (v, 0). 

Now observe that {((/)"(u), — n)} describes both the random walk and 
the tree starting from (u, 0). Also if A^^ n A^^ = 0, then ttIv = Iv and, 
more importantly, C(v) = ip{v)- Hence the "random walk" and the tree from 
(u,0) are coupled and so are the "random walk" and the tree from (v,0). In 
particular, this happens when both < [||u — v||i/2] and my < [||u — v||i/2]. 
Let ko = ||u — v||i/2. From the above discussion, we have 



for constants Ci = 2 and C2 = (l/2)|log(l —p)\- 

With the above estimate at hand, we look at the process {(0"(u), C"(v)) : n > 
0}. Without loss of generality we take u = o. For e > and constant K > 
(to be specified later), define 



lP({C(v) / V(v)}) < P({(C/i"(z)) >p for all z G (u + A^J} 

U {{UU^)) > P for ah z G (v + A^J}) 
= 2P({(C/{'(z)) > p for all z G (u + A^J}) 
= 2(l-p)#'^'=o. 

Since (1/2)A;^ < t^A^ < 2k^, the above inequality gives 
(17) P({C(v) = V'(v)}) > 1 - Ci exp(-C2||u - v||?) 



(18) 



S„,,(v) := {C'(v) G <A"'(o) + (A„2(i+.) \ A„2(i-.,), 

IIC(v) - ^\o)\\i > Klogn for alH = 1, . . . ,n^}. 
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This event is an independent random walk version of the event An,s{v,0) 
defined in (14), except that here we require that the two random walks come 
no closer than Klogn at any stage. 

We will show that there exists a > such that 



(19) 



sup 

ve(A,^(i+s)\A„(i-.)) 



,(v))^)<C3n-" 



for some constant C3 > 0. 

Since (5„,s(v))^ C ^„,^.(v) U F„,^.(v) U Gn,e{y), where 

EnA^) := {IIC(v) - cP\o)\\i < Klogn for some i = 1, . . .,71"}, 

to prove (19) it suffices to show the following. 

Lemma 3.3. There exist a > and constants 0^,0^,0^ > such that, 
for all n sufficiently large, we have: 

(a) sup,g(A^^j^^^,\A„^,_,))IP(^n,.(v)) < Qn-", 

(b) sup,,(^^^,^^^,\^^^^_^,)P(F„,,(v)) <C5n-, 

(c) sup,,(^^^,^^^,\^^^^_^,)P(G.,,(v)) <C76n-^ 

Proof. First we fix v G (A„(i+,) \A„(i-,)). Since {((/>'' (o), C(v)) : n > 0} 
and {(X]r=i ^i^^ + Z]r=i ^i) '■'n> 0} have the same distribution, we have 



< K log n for some i = 1, . . . , n 



E ~ E ^ (v + Ai^iogn) for some i = 1, 



,71 



<^\Y.^3-Y. ^ + logn) for some i>l\ 

= P [ U I ^ - ^ = z for some « > 1 M • 

Now ^*-^^(Xj — Yj) is an aperiodic, isotropic, symmetric random walk 
whose steps are i.i.d. with each step having the same distribution as X — y, 
where Y is an independent copy of X. Since Var(X — Y) = 2 Var(X) = 2a^I 
[where cr^ := Var(X(l)) and Var(X) denotes the variance-covariance matrix 
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of X] and Eugz3 |upP(X - y = u) < oo, by Proposition P26.1 of Spitzer 
[(1964), page 308], 

(20) lim |z|p| VXj -VFj =z for some i> ll = (47rVar(X(l)))"\ 

For V G (A„{i+e) \ A„{i-£)) and z E v + Aft-iogm we must have that, for 
all n sufficiently large, |z| > n^~^/2. Thus for all n sufficiently large and for 
some constants Ci,C^^C^ > 0, we have, using (20), 

F{En,e{y)) < J2 ^\j2^3-J2^j=^ foi' some i > 1 1 
ze(v+AKiog„) li=i j=i ) 

This completes the proof of Lemma 3.3(a). 

For the next part of the lemma, observe that, for sufficiently large n and 

all vG A„(i+,) \A„(i_,), 



4 

2{l+e) 



(21) 



> n 



< 



< 



>„2(l+.)_^(l+.) 



To estimate the above probability let X — y = Z = (Z(l), Z(2), Z(3)), where 
E{Z{i)) = and Var(Z(i)) = 2a^. Then, letting Ei=i(^i - denote 
the ith co-ordinate of the process J2^=i{^j ~ ^j) ^i^d using Chebyshev's 
inequality, we have 



> 



n 



2{l+e) 



< 



.i=l 



T.(^j-yj)i^ 



> 



n 



2(l+£) 



6 



< 3P 
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> 



6 



< 



< 



3n''Var(Z(l)) 

(„2{l+£)/Q)2 

Cii 



n 



Ae 



for some constant Cn > 0. Combining the above inequality with that ob- 
tained in (21), we have 



sup W{Fn,e{v))<^, 

ve(A^(i+,)\A^(i_,)) 



which proves Lemma 3.3(b). 

Finally, for the last part of the lemma, we have that if 0<e< 1/3 and 
V G A„{i+£) \ A^{i_e), for all sufficiently large n, ||v||i < n'^^^~'^\ Therefore, 



"(Gn,.(v))< 



(22) 



< 



< 



< n 



2{l-e) 



< ||v||i + n 



< 



2n2(i--) 



i=l 



< 



2^2(1-.) 



By the central limit theorem, as n ^ oo, ^^=i{Xj — Yj){l) / (\/2an'^) con- 
verges in distribution to a random variable N (say) with a standard normal 
distribution. Thus 



(23) 



< 



E^=i(x,-y,)(i)| ^ V2n-2^' 



a/2, 



< 



3cr 



V2n-2^ 



3o- 

|iV|< 



V2n- 



2e 



3cr 
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Of the terms in the above inequahty, we have 



(24) 



3a 



< 



V2n-2^(3o-)-i V2vr 



exp i ——] dx 



2tt 



and we use Berry-Essen bounds [see Chow and Teicher (1978), Corollary 
9.4, page 300], to obtain 



< 



3(T 



P \N\ < 



2e 



3(T 



< 



< 



3cj 



N< 



2e 



(25) 



+ 



<2supP<i 



3a J 

<x \ -P{N<x} 



3a 

N<- 



■2e 



3a 



< 



Ci2E{Zf) 



for some constant C12 > 0. Combining (22)-(24), we have Lemma 3.3(c). □ 

Proof of Lemma 3.2. Let v := (v,0) G Z^. Observe that An,e(v) 5 
BnA^) n {5^0) = Ei=i Xj, g'{v) = V + E;=i Yj for all 1 < i < n^}. Hence 



> P B„,,(v) n 5^(0) = ^ Xj, g\w) = V + ^ for 1 < i < 



n 



i 



Sn,.(v) n g*(0) = 5^ Xj, <7^(v) = V + 5^ for 1 < i < - 1 
I j=i j=i 



X P 5"' (0) = E ' fl"' (v) = V + E 



Sn,.(v) n (7'(0) = ^ X,, g*(v) = V + ^ y,- for 1 < i < - 1 
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> p| B„,,(v) n L^(0) = ^ X.j,g\^) = V + ^ for 1 < ^ < - 1 
I I i=l j=l 

X (l-Ciexp(-C72(/^logn)3)), 

where the last inequahty follows from (17) after noting that, given i3„,e(v), 
gi[Q) = Ej=i^j and 5*(v) = V + Ej=i^i hold for all 1 < i < - 1, we 
have ||5f"^~^(0) — fi'"'''^^(v)||i > J^logn. The above argument may be used 
iteratively for i = 1, . . . , — 1, and together with (19), we have 

P(A,e(v)) > (1 - Ciexp(-C2(i^logn)3))"V(S„,,(v)) 

> (1 - Cin^exp(-C2K^logn))(l - Cgn"") 

> (1 - Cin^n-^^^''){1 - Can-") 
= (1 - Cin-^2i^'+^)(l - Cgn""). 

Taking K such that C2K^ > 4 [i.e., > 8| log(l — p)\ we have 
P(A,e(v)) > 1 - Cin-^2^'+4 - Can-" 
> 1 - Cn-^, 

for some constant C > and /3 := min{a, C2K^ — 4} > 0. This completes the 
proof of Lemma 3.2. □ 

Finally, to complete the theorem we need to show that Q admits infinitely 
many trees almost surely. For k>2, define D^{n, e) := {(ui, U2, . . . , u^) : Uj e 
such that n^~'^ < ||g*^(uj) — (7''(uj)||i < n^^^ for all i ^ j}. Define the event 
A(n,e,Ui,U2,...,Ufc) :=K(i--) < \\g^\ui) - g^\u,)\\i < n^(^+'^ andg*(u,) / 
g*(uj) for alH = 1, . . . , and for all i / j}. Using Lemma 3.2, we can easily 
show, for < e < 1/3 and for all large n, 

(26) inf{P(A(n, e, ui, U2, . . . , Ufe) : (ui, U2, . . . , Ufc) G L>^-(^, £))} > 1 " 

where is a constant independent of n (depending on k) and /3 is as in 
Lemma 3.2. We may now imitate the method following the statement of 
Lemma 3.2 to obtain 

P{/(ui) / g\uj) for ah t > 1 and for 1 < i / j < /c} > 0. 

Thus, by translation invariance and ergodicity, we have that, for all k>2, 

¥{Q contains at least k trees} = 1. 

This shows that G contains infinitely many trees almost surely. 
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4. Geometry of the graph Q. We now prove Theorem 2.2 for d = 2; with 
minor modifications the same argument carries through for any dimensions. 
The idea behind this proof was suggested by the referee. 

For t £ Z, consider the set Nt := G Ci {y = t}, the set of open vertices 
on the fine {y = t}. For x £ Nt and n > 0, let B^{x) := {y £ Nt+n ■ h"-{y) = 
x} be the set of the nth-order ancestors of the vertex x £ Nf. Now con- 
sider the set of vertices in A'^ which have nth-order ancestors, that is, 
Mj^"^ := {x£Nt: 5f (x) / 0}. Clearly, iV^"^ C iV^"^ for n > m and so Rt : = 

(n) (n) 

lim.„^oo Aft = r\n>o^H is well defined. Moreover, this is the set of ver- 
tices in Nf which have bi-infinite paths. We want to show that ¥{Rt = 0) = 1 
for all t£TL. Since : i G Z} is stationary, it suffices to show that P(-Ro = ^) 
= 1. 

First note that by the translation invariance of the model, P({7^i?o = 
0} U {#-Ro = oo}) = 1- Now suppose P(#-Ro = oo) > 0. A vertex x £ Rt is 
called a branching point if (x) n Rt+i) > 2, that is, x has at least two 
distinct infinite branches of ancestors. Note that this notion of "branching 
point" is similar to that of "encounter point" of Burton and Keane (1989). 
As in their proof of the uniqueness of the percolation cluster, our proof 
essentially uses the fact that it is impossible to embed a tree in a lattice. 

We first show that 

(27) P(Origin is a branching point) > 0. 

Since F{#Rq = oo) > 0, we may fix two vertices x = (xi, 1) and y = {yi, 1) 
such that 

¥{x,y£iBliO)nRi))>0. 

Thus the event Ei := {B^{x) ^ 0,B'^[y) ^ for all n > 1} has positive 
probability. Further, this event depends only on sites {u := (ni,U2) :U2 > 1}. 
Now, consider the event £2'.= { (i, 0) is closed for alH 7^ with —2 max{ |xi | + 
l,|yi| + l} < i < 2max{|2;i| + l,|yi| + l} and (0,0) is open}. Clearly P(£;2) > 
0. Since Ei and E2 depend on disjoint sets of vertices, we have 

P(Origin is a branching point) > P(£;i n E2) = F{Ei)¥{E2) > 0. 

Now, we define ro(n) := #(-Ro H ([— n,n] x {0})) and ri(n) := #{Ri Pi 
([— n, n] X {!})). We arrange the points of i?on ([— n, n] x {0}) as ui, . . . , tiro{n); 
in an increasing order of the x coordinates. By our construction of Q, neither 
U2 nor n^(,(„)_i nor any of the vertices between them can be connected to 
a vertex on A^i which lies outside [— n, n] x {!}. Thus, each of the vertices 
U2,U3, . . . ,u.ro{n)-i '^iH have at least one ancestor in the set Ri H ([— n,n] x 
{!}). Moreover, each of the branching points in U2, . . . ,UrQ{n)-i has at least 

two distinct ancestors in the set i?i n ([— n,?i] x {!}). Thus, if Tq '{ii) is the 
number of branching points in [— n,n] x {!}, we must have 

(28) ri(n)-(ro(n)-2)>rJ')(n)-2. 
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But, by stationarity, we have E(ri(n)) = E(ro(n)) for all n > 1. Thus, for n 
sufficiently large, from (27) we have 

= E(ri(n) -ro(n)) >ErJ^^(n) -4 
= {2n + l)P(Origin is a branching point) — 4 > 0. 
This contradiction establishes Theorem 2.2. 



5. Limit theorem. We first prove Theorem 2.3(a). The proof of the next 
part of the theorem is similar and thus omitted. For simplicity in notation 
we shall prove the result for d = 2; however, our method is also valid for 
higher dimensions. 

Fix 1/ >0. Let Bn := [i-,n] x [l,n] be a box of width n and, for (i, j) € 
Bn n I?, define random variables Yij as 

Y f 1, if the degree of the vertex (i, j) in i?„ n V is + 1, 

^'-^ ' \ 0, otherwise. 

Note for a vertex (i, j), Yij = 1 if and only if there are exactly v edges "going 
up" from and one edge "going down" from it. 

Let y/"^ := Er=i(^i,i - HYid)) and Sn := Ei=i ■ To prove Theorem 
2.3 we need to show that the distribution of Sn/n is asymptotically normal. 

Towards this end, first observe that, for fixed j, {i^j}i>i is an a-mixing 
sequence of random variables; that is, for all m > 1, A G cr(Yi j , Y2,j, . . . , Ymj) 
and B e a{Y^+nj,Ym+n+i,j, ■ ■ ■), we have \F{AnB) -¥{A)F{B)\ < a„, where 
a„ — > as n — > oo. Indeed, given A and B as above, define 



E := <^ there exists an open vertex in each of the sets 

f , , n 3n] ( 3n n 1 

Ui,j):m + j<i<m + — '>,Ui,j + l):m + — <i<m+ -y 



f n 5n] ( 

|(z,j + l):m + -<i<m+ — |,|( 



3n 3n 
+ — <i<m + — 



Now ¥{E) = (1 - (1 - p)"/8)4 ^ 1 as n ^ oo. Also, given E, the event A 
depends only on the configuration of the vertices — < m + ^}, 

{(hj) - "i- ^ "i} and + 1) :i < m + ^}, while the event B depends on the 
vertices — 1) : i > m + ^}, :i>m + n} and + 1) :% > m+ ^}. 

These sets of vertices being disjoint, given E, A and B are conditionally 
independent, a simple conditioning argument now yields that, for n large 
enough, 

\F{A r\B)- F{A)F{B)\ < SFiE") < Ci exp(-C2n) 

(29) 

for constants Ci, C2 > 0. 
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Also observe that, for fixed i, {^ij}j>i is a one-dependent sequence of 
random variables; that is, for fixed i, Yi j is independent of l^j' for j' ^ 

3 - 1, j,j + 1- 

Now, for some < 6 < 1 to be chosen later and for < k < Vn, where 

._yW J L y(") 

(n) 

-'^(fc+l)[n'5J+fc+l' 

First we show that, for any r > 1, there exists a constant C > such that 

(30) E(y/") + • • • + )^ < Cr2?i2. 

Indeed note that, as in the proof of the first part of Theorem 27.5 of 
Billingsley (1979), we have E(y/"^)^ = E(y/"^)^ < Kn^ for some constant 
K>0. Now 

(31) E(j2"^lA = E E(i"i"^y;^"Vwy/")), 

\fc=l / k,l,s,t=l 

and using the fact that {^fc"^}fc>i is a one-dependent sequence of random 

~ (n) 

variables, the Cauchy-Schwarz inequality and that =0, we obtain 

after some elementary calculations 



\fc=i / 

Here the term 2rE(Y'/"^)^ comes from the terms in the sum J2'j,k,s,t=i^0^j"^^'^k'^'^ ^ 
yi"V/"^) when j,k,s,t are close to each other so as to have dependence 
among all the four random variables making the product, while the term 
comes from the terms of the sum when j, k are close to each other, 
s,t are close to each other, but there is independence between (l^'-"^ Y"^^"^ ) 
and (yi"\y/"^). This proves (30). 

Now taking r= [n^\, and using the fact that Wg^"^ . . . are i.i.d. 

random variables, we have from (30) that E(W^''y < Cn2+2^ for all /c > 1. 
Also 



/ J \ 
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(32) =[7i'\E{Yl''^f + 2 J2 Cov(y/"\y/;^i) 

i=i 

= [n^jE(y/"^)V2([n^J - 1) Cov(y/"\ Fa^"^). 
In the above expression, 

n—l n—s 

E(y/"V = nVar(yi,i) + 2 J] J] Cov(y,,i, y,+t,i) 

s=l t=l 
n-1 

= n Var(yi,i) + 2 ^(n - s) Cov(yi,i, yi+,,i) 

s=l 

= 0{n) as n — > cx), 

where the last equahty follows because from the a-mixing of the sequence 
{Yt,i}t>i we have X]t^2 Cov(yi^i, y^^i) < CY^'^2'^t < oo for some constant 
C > 0. Moreover, by the Cauchy-Schwarz inequality, 

Cov(y/"\y2("))<E(y/"V- 

Thus, from (32), we have Var(VF^^"'^) = 0{n^^^) as n ^ oo and 

=Drr,(i-'5)+(i+'5) 

Finally, for < 5 < 1, 



(33) Var =0(n(i-'^)+(i+'5)) =0(n2) asn^oo. 



lim V ^E(W^("))^ 

™,ei(VarE^=iW^f¥ 

- 2+25 



< lim y C — -r- = lim Cn"^"! = 0. 

fc=l 

Thus by Lyapunov's central limit theorem [see Billingsley (1979), Theorem 

27.3, page 312] we have that, for < 5 < 1, l/(\/Efc=i Var(W^"^)) Y.l=i W^i"^ 
converges in probability to a standard normal random variable. 
Now let 7]^ := Z]fc=i ^i"^- We will show that 

(34) rjn/n^Q in probability as n ^ oo. 

Indeed, 



E(4"V < EVar(yi,fc) + 2n^Cov(yi,fc,y,,, 
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oo 

< ?iVar(yi,i) + 2n^Ciexp(-C2i) 

i=2 

< Mn for some constant M > 0. 
Thus, using the fact that r„ = 0{n^~^) as n — > oo, we have, for e > 0, 

;> nt) :^ — 
This proves (34). 

To complete the proof, we have to show that — > in probabihty as 
n — > oo. First observe that number of terms in En is at most [n^J . Therefore 
taking 5=1/2, from (30) we have E(£'^) < Cin?. Hence, for e > 0, 

(35) P(|£;„| >ne) < -^^^0 asn^oo. 

Theorem 2.3(a) now fohows by combining equations (34) and (35) and the 
fact that X]fc=i ^A:''"V'^ asymptotically a A^(0, s^) distribution, where 



?7„|>ne)< — = ^0 asn^oo. 



= Var(yi,i) + 2 ^ Cov(yi,i, 



i=2 



Y 



+ 2^ Cov(yi,i, y,,2) + 2 £ Cov(yi,2, y^,i). 

i=\ i=2 

Note that to compute we use the fact that {(y^ j, yjj_|_i)}j>i is an a- 
mixing sequence. 

6. Degree of a vertex. To prove Proposition 2.1, observe that, given the 
vertex (0,-1) is open, let 

1, if the vertex (0,0) is open, 
0, otherwise, 

Xi = 0) : i < — 1 : (i, 0) is connected by an edge to (0, —1)}, 

X2 = 0) ■i> 1 : («, 0) is connected by an edge to (0, —1)}. 

Clearly the degree of (0,-1) equals Y + Xi + X2. Now given the vertex 
(0, —1) is open, the probability that the vertex (— /, 0) is connected to (0, —1) 
and that there are exactly r — 1 vertices in {(i, 0) : — / + 1 < i < — 1} which 

are connected to (0,-1) equals (^^_^^p'^(l — p)'~''(l — p)^'~^((l — p) + ^p)- 

Thus F{Xi >r)= {lz\)p'-{l-py-'{l-pf-H{^ -P) + hp)- An easy 
calculation now completes the proof of the proposition. 

Similarly, in two dimensions, given that a vertex v is open, the distribution 
of the number of edges of length / "going up" from v is binomial with 
parameters 2 and (1 — |)(1 — p)^'~^. 
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Remark 6.1. From the above distributions we may calculate the quan- 
tities K{Sn), Var(5n), and the related quantities involving L„ required in 
Theorem 2.3 for two dimensions. 

Acknowledgments. We are grateful to Professor S. Popov for his sugges- 
tions regarding the proof for d = 3 and to an anonymous referee for suggest- 
ing that the Burton-Keane argument would yield Theorem 2.2. 
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